Wilsonian renormalization group analysis of nonrelativistic three-body systems 

without introducing dimerons 
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Low-energy effective field theory describing a nonrelativistic three-body system is analyzed in 
the Wilsonian renormalization group (RG) method. No effective auxiliary field (dimeron) that 
corresponds to two-body propagation is introduced. 

The Efimov effect is expected in the case of an infinite two-body scattering length, and is believed 
to be related to the limit cycle behavior in the three-body renormalization group equations (RGEs). 
If the one-loop property of the RGEs for the nonrelativistic system without the dimeron field, which 
is essential in deriving RGEs in the two-body sector, persists in the three-body sector, it appears to 
prevent the emergence of limit cycle behavior. We explain how the multi-loop diagrams contribute 
in the three-body sector without contradicting the one-loop property of the RGEs, and derive the 
correct RGEs, which lead to the limit cycle behavior. The Efimov parameter, so, is obtained 
within a few percent error in the leading orders. We also remark on the correct use of the dimeron 
formulation. 

We find rich RG-flow structure in the three-body sector. In particular, a novel nontrivial fixed 
point of the three-body couplings is found when the two-body interactions are absent. We also find, 
on the two-body nontrivial fixed point, the limit cycle is realized as a loop of finite size in the space 
of three-body coupling constants when terms with derivatives are included. 



I. INTRODUCTION 

Nonrelativistic three-body systems have many interest- 
ing features. There have been long-standing discrepan- 
cies and anomalies between theoretical calculations and 
the experiments for several physical quantities in the 
three nucleoli systems. (See Ref. [l| for a review.) So- 
called Borromean nuclei, three-body bound states, any 
two of whose constituents do not form a bound state, 
arc known to exist. (See Ref. @ for a review.) V. Efi- 
mov showed that there are infinitely many three-body 
bound states with the ratios of the binding energies of 
the subsequent bound states being a universal constant, 
E n+l /E n = e - 27r / s ° = 1/515.03 • • • , when the two-body 
scattering length is infinite 0, E| ■ This "Efimov effect" 
has recently been attracting much attention because of 
the advancing experimental control over the very cold 
atomic systems by using Fcshbach resonance. 

At low energies, the internal degrees of freedom of the 
particles (i.e., nucleons, nuclei and atoms) become irrele- 
vant and the system is well described by an effective field 
theory (EFT), in which only the particles without inter- 
nal structure interact locally. EFT is an efficient, model- 
independent approach and admits a systematic improve- 
ment of the description. 

Bcdaque, Hammer, and van Kolck Q first noticed that 
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the Efimov effect is related to the RG limit cycle behav- 
ior of the three-body nonderivative coupling of the EFT. 
RG analysis has played an important role in recognizing 
that the Efimov effect is a new kind of universal phe- 
nomena in the three-body systems [(| . There are several 
papers which deal with the Efimov effect in the light of 
nonperturbative RG analysis. For a recent review, see 
Ref. 0. 

In their analysis, as well as subsequent studies by other 
authors, the so-called dimeron field Q, an auxiliary ef- 
fective field which represents two-body propagation, is 
exploited. To our best knowledge, there is no literature 
in which the connection between the limit cycle behavior 
and the Efimov effect is established without a dimeron 
field. One might wonder if the limit cycle behavior can 
be obtained without introducing a dimeron field. 

Although dimeron is a useful device, its introduction 
needs some care. It is difficult to establish the relations 
between a set of the coupling constants in the theory 
with dimeron and that in the original theory because of 
a nontrivial factor in the path-integral measure arising 
from the field transformation, an awkward contribution 
which depends the regularization of how to define a func- 
tional determinant. (In Appendix lAl we explain how the 
additional contribution from the Jacobian affects the re- 
lations among the couplings.) It is thus easier to think 
the theory with dimeron as another EFT, the couplings 
of which are to be determined by a matching procedure, 
than to make a direct connection between the two theo- 
ries. 

In this paper, we perform an RG analysis for a nonrel- 
ativistic three-body system of a single (complex) scalar 
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field without introducing the dimeron. (It is straightfor- 
ward to extend our analysis to the fcrmionic case.) The 
RGEs are derived and the limit cycle behavior of the 
three-body coupling is identified for the first time in this 
formulation. We also find that the usual formulation with 
dimeron field misses a certain contribution, which seems 
essential for obtaining the Efimov parameter sq within a 
few percent error. 

It is well-known that nonperturbative RGE has the 
"one- loop" property, i.e., the running of the coupling con- 
stants is determined by a sum of one-loop diagrams with 
the propagator being the full one. In a nonrclativistic 
system, because of the particle number conservation and 
the absence of the antiparticle, there is no dressing for 
the propagator. This nonrelativistic feature divides the 
space of states into sectors labeled by particle numbers, 
and makes the RGEs in the two-body sector very simple: 
one only needs to consider actual one-loop diagrams that 
satisfy the particle number conservation, with the free 
propagator. 

If this one-loop property (with free propagator) of the 
RGEs persists in the three-body sector, there appears to 
be no chance for the RGEs to have a limit cycle solution. 
This is a puzzling situation: on one hand, the Efimov ef- 
fect implies the limit cycle behavior when the two-body 
couplings are tuned to the critical values, on the other 
hand, the general consideration given above implies the 
one-loop property, which seems to prevent the limit cycle 
behavior. It is one of our main results to explain how the 
one-loop property of the RGEs is consistent with con- 
tributions from multi-loop diagrams in the three-body 
sector. The resulting RGEs do allow limit cycle behav- 
ior. 

We find very rich structures in the RG flow of the 
three-body coupling constants. In the leading-order cal- 
culation, in which only the nonderivative three-body cou- 
pling is included, we find a nontrivial fixed point on the 
two-body trivial fixed point. It implies the existence of a 
three-body strong coupling phase with the phase bound- 
ary on which the nontrivial fixed point resides. 

Extending the space of three-body operators up to in- 
cluding the terms with two derivatives, we find that, on 
the two-body nontrivial fixed point, the limit cycle is re- 
alized as a loop of finite size, and a nontrivial fixed point 
associated to it. 

We also emphasize that the nonzero values of the two- 
body effective range and the off-shell parameter, which 
does not affect the on-shcll two-body amplitude, mod- 
ify the the ratios of the binding energies of the three- 
body bound states. Although the effects of the nonzero 
value of the two-body effective range have been investi- 
gated [H , those of the off-shell parameter have never 
been considered. 

Several comments on the literature are in order. Al- 
though Bedaque, Hammer, and van Kolck Q first showed 
the connection between the RG limit cycle and the Efi- 
mov effect, it does not come as a direct consequence of 
the RGE analysis, but rather from the scale-invariant 



behavior of the three-body amplitude. Actually they 
obtained the cutoff dependence of the coupling of the 
three-body nonderivative contact interaction after de- 
termining the Efimov parameter. Moroz, Floerchinger, 
Schmidt, and Wetterich ll| derived the Efimov param- 
eter in the functional RG formulation with the dimeron 
field. Their leading-order value of the Efimov parameter 
contains about 40% errors. Their numerical calculation 
for the full RGE shows, however, that the resulting value 
of the Efimov parameter is in good agreement with the 
Efimov's. 

There are several other papers devoted to the RG anal- 
ysis of the three-body systems. For example, Dichl, 
Krahl and Schrerer [12| study the ratio of dimeron- 
particle scattering length of to that of particle-particle. 
Krippa, Boris, Walet, Niels and Birse [13| do a similar 
thing for the ratio of dimeron-dimeron scattering length 
to that of particle-particle. All of these employ dimerons 
in their analysis. 

The paper is organized as follows. In Sec. UH we re- 
capitulate the importance of Wilsonian RG analysis in 
revealing the physical features of few-body systems de- 
scribed by EFT. The one-loop property "puzzle" is ex- 
plained and the solution is given in Sec. IIIII Wc derive 
the RGEs for the coupling constant for the three-body in- 
teraction in the leading order in derivative expansion in 
Sec. lIVI Comments on the difference between the dimeron 
formulation and ours are also given. In Sec. |Vj we first 
concentrate on the case in which the two-body effective 
range and the off-shell parameter are zero. The RG flows 
are given and the Efimov parameter is obtained on the 
two-body nontrivial fixed point. The other cases are also 
examined. A novel nontrivial fixed point in the three- 
body coupling is identified. We examined the next-to- 
leading order corrections in Sec. I VII The RG flow in the 
three-body coupling space is presented. The limit cycle 
is found to be realized as a loop of finite size in this space. 
We summarize the results in Sec. IVIII 

In Appendix [A] we illustrate a systematic way of in- 
troducing the dimeron field in the path integral formu- 
lation. We argue that there is a nontrivial Jacobian in 
the measure due to the nonlinear change of variables. It 
leads to change of coupling constants which has been ne- 
glected in the literature. In Appcndix[Bj we demonstrate 
the derivation of the RGEs by using an example. In Ap- 
pendix [D] we explain that the pole appearing in the RGE 
for the three-body coupling is due to the existence of the 
two-body bound state. 



II. EFFECTIVE FIELD THEORY IN THE 
THREE-BODY SECTOR AND WILSONIAN RG 



In this section, we recapitulate the basic idea of 
EFT [3 and the usefulness of Wilsonian RG [lj] in the 
context of EFT. 
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EFT 



the naive dimensional analysis. 



Field theoretical description of nature has a hierarchi- 
cal structure: at every scale there are relevant degrees of 
freedom and the interactions among them are described 
in terms of local operators subject to the symmetry of the 
system. Even though the system is composed of com- 
posite particles, if the momentum scale in question is 
smaller than the scale of the internal structure, one may 
neglect the structure and describe the system in terms 
of the fields corresponding to the particles. The effects 
of heavier particles that are not included in the EFT are 
encoded in the values of the coupling constants of local 
operators. Such a field theory is called an effective filed 
theory (EFT). A physical system may be described by 
various EFTs at various momentum scales. 

The EFT description is very general: an EFT is defined 
by the relevant degrees of freedom, dimensions of space- 
time, and the symmetries. It contains all possible opera- 
tors with these attributes, which are infinitely many. A 
single EFT may describe several different systems that 
share the same attributes. Features of a specific system 
are reflected in the values of the coupling constants. Thus 
the EFT description is model independent. 

The accuracy of the EFT description is controlled by 
power counting rules. Power counting rules tell us the de- 
grees of importance ("orders") of operators. The higher 
the order is, the less important the contributions are. To 
a given order, one only needs to consider a finite number 
of operators and a restricted set of diagrams. If the power 
counting rules are consistent, counterterms are also sup- 
plied to the order. Thus the EFT description is system- 
atically improvable. 



B. Wilsonian RG reveals nonperturbative aspects 
of EFT 

In several cases, nonperturbative physics shows up in 
the EFT context. It is of central importance to establish 
the power counting rules for such a case to extract phys- 
ical information out of EFT. An example is the nuclear 
effective field theory (NEFT) [H-G]| for the two-nucleon 
system in the S waves. The existence of a bound state, 
deuteron, in the spin-triplet channel, is a clear sign of 
nonperturbative dynamics. 

It is the Wilsonian, or nonperturbative, RG analy- 
sis [HI, [l9l - l22j that reveals properties of the nonperturba- 
tive dynamics. The effects of quantum fluctuations are 
examined scale by scale by changing the cutoff in the 
Wilsonian RG analysis, without relying on the perturba- 
tive approximation. 

The application of such analyses to the NEFT in two- 
nucleon sector is given in Refs. (23l - [25j . The nonpertur- 
bative feature of the system is translated into the RG lan- 
guage as the existence of a nontrivial fixed point. Because 
of the large anomalous dimensions, the power counting 
rules arc modified from the pcrturbative ones based on 



C. Efimov effect in the light of RG 

More than forty years ago, V. Efimov considered the 
infinite two-body scattering length limit of a nonrelativis- 
tic system, probably motivated by the large two-nucleon 
scattering lengths in the S waves. He started with the 
Schrodinger equation with a short-range two-body po- 
tential and derived an effective three-body Schrodinger 
equation. He then noticed that there are an infinite num- 
ber of three-body bound states with the ratios of binding 
energies of subsequent states being a universal constant. 
The result does not depend on the details of the short- 
range two-body potential. 

One may think that this phenomenon can be explained 
in the language of RG in an EFT. Such a formulation 
would make the universal feature of the phenomenon 
more transparent. Let us consider the EFT Lagrangian 
of a single nonrelativistic boson of mass M: 
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(2.1) 



where and the ellipses stand for 

the terms with more derivatives. Throughout this paper, 
we concentrate on the S waves, so that the operators 
contributing to higher partial waves are not included. We 
do not consider the operators which act on more than 
three-body states, e.g., ((^) 4 </> 4 > (</^) 5 5 i etc., because 
we are interested in the three-body sector and they do 
not contribute to the sector. Note that we have included 
the so-called redundant operators such as the one on the 
third line in Eq. (|2.1J). w hich are necessary to renormalize 
the theory off shell [26|. 

In an important paper, Bedaque, Hammer, and van 
Kolck Q use another version of EFT with the so-called 
dimeron field, D. In the leading order the Lagrangian is 
given by 
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They find the limit cycle behavior for g% from the scale 
invariant property of the three-body amplitude in the 



limit of an infinite scattering length, 
A 2 52 (A) 



H(A) 



4M 5 ?(A) 
sin [so ln(A/A„) — arctan(s ( j" 1 )] 
sin \sq ln(A/A*) + arctan(sg 1 )1 



(2.3) 



where sq is called the Efimov parameter, sq = 
1.00624 • • • , A is the floating cutoff, and A* is a constant. 

Note, however, that although they find the intimate 
connection between the Efimov effect (as a consequence 
of discrete scale invariance) and the RG limit cycle be- 
havior, they obtain the (floating) cutoff dependence of 
the coupling constant not from the RGEs. Note also 
that their results are given only for the case of an infinite 
scattering length. 

More recently a field-theoretical derivation of the limit 
cycle behavior directly from the RGEs is g iven by Moroz, 
Flocrchinger, Schmidt, and Wetterich [Tl[. They also use 
the dimeron field. Their formulation admits an arbitrary 
momentum dependence of the coefficient functions. In 
the leading-order approximation (the point-like approx- 
imation), the RGEs are solved analytically, but in the 
higher orders they are solved only numerically. 

Even though the use of dimeron is useful, it is just an 
option. The RGEs and the limit cycle behavior should 
be obtained even if the dimeron field is not introduced. 
Such a demonstration would provide a better insight into 
the use of dimeron field. 



III. THE ONE-LOOP PROPERTY PUZZLE 

In this sections, we explain the one-loop puzzle and 
its solution. For simplicity, we consider a nonrclativistic 
bosonic system described by the Lagrangian (|2.1[) . The 
dimeron field is not introduced. 



A. one- loop property of nonperturbative RGEs 

It is well known that nonperturbative RGEs have a 
kind of "one-loop" property. The Wegner-Houghton 
equation [l9| for the Wilson action S e g, for example, 
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(3.1) 



has the "one-loop" term (the first line) , together with the 
"dumbbell" term (the second line) , where the integration 
with prime means the integration over shell mode, 1 — 
dt < p < 1 . In the functional flow equation [2l|, [22j , 
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FIG. 1: Three types of time-ordered one-loop diagrams which 
are apparently the only diagrams contributing to the running 
of three-body coupling constants. 



where Tk[4>] is an averaged action with the momentum- 
shell parameter k, on which the regularization function 

(2) 

Rk depends, T k is the second derivative of Tk with re- 
spect to 4>, and t = ln(/c/Ao), the whole contributions 
come from the "one-loop" diagrams. In both cases, the 
"one-loop" diagrams are composed of the full propaga- 
tors, so that the actual structure is much more compli- 
cated than it appears. 

In nonrelativistic systems, because of the absence of 
antiparticles, the "one-loop" diagrams are really one- 
loop, i.e., the propagators are bare ones. There are no 
tadpole-type diagrams, in which an internal line starts 
and ends at the same vertex. 

In the two-body sector, there is only one type of one- 
loop diagrams, with various kinds of vertices. It gives 
rise to the RGEs of the couplings of two-body contact 
interactions. Note that, since the two-body amplitude is 
given as the sum of bubble chains, the cutoff indepen- 
dence of the two-body amplitude also leads to the same 
RGEs. 

For the theory under consideration, if the one-loop 
property persists in the three-body sector, there appear 
to be only three types of one-loop diagrams that con- 
tribute to the running of three-body coupling constants. 
Sec Fig. [TJ It is clear that the RGE for do would not 
exhibit periodic behavior. To be concrete, the RGE for 
do would be 



dv 

— = -A.v - 
dt 



ax 



(3.3) 



where x and v are dimensionless coupling constants de- 
fined by 



MA 
4^2 



Co, 



MA 4 
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(3.4) 



and t = ln(Ao/A) with A being the floating cutoff, a 
and (3 are positive dimensionless constants. In the limit 
of infinite scattering length, x —¥ — 1, the RGE (|3.3[) 
can be easily solved but the periodic behavior cannot 
be obtained. It seems that multi-loop diagrams must 
contribute to the RGE for do. But how? 



B. normal-ordering 

It is important to note that, in nonrelativistic systems, 
the interaction operators are all (implicitly) assumed to 
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FIG. 2: A time-ordered diagram which contributes to three- 
body local counterterms when the loop momentum is set 
in the shell. The resulting counterterms are already in the 
normal-ordered form. 



be normal-ordered. Under this assumption, we can con- 
sider the sectors labeled by the particle number sepa- 
rately, and the n-body operators, such as ($) n (j) n do not 
affect the m-body sectors with m < n. This is a very 
favorable feature we would like to keep to simplify the 
calculations. 

The crucial point is that the counterterms generated by 
RG transformations are in general not normal-ordered. 
In order to obtain the correct RGEs for the coupling 
constants for the normal-ordered operators, one needs 
to rewrite the counterterms in the normal-ordered form. 
This rewriting of an n-body couterterm leads to normal- 
ordered m-body operators with m < n. 

Let us explain what is going on by examples. The first 
example is the diagram shown in Fig. [2] With the loop 
momentum being in the shell A— 5 A < p < A, it generates 
an effective local interaction which can be canceled by 
3-body counterterms that are all normal-ordered. Note 
that the diagram is a time-ordered one. It is impossi- 
ble to contract any external outgoing and incoming lines 
without going backward in time. If we assign an exter- 
nal outgoing line and an incoming line with the creation 
and annihilation operators, and a, respectively, this 
diagram has the (a ! ) 2 (a' ! a)a 2 structure, which is already 
in the normal-ordered form. On the other hand, in the 
second example shown in Fig. [3] with the loop momen- 
tum being in the shell, the diagram generates the con- 
tributions that are not normal-ordered. In this case, it 
is possible to contract the external lines without going 
backward in time. In other words, the diagram has the 
(a') 2 '(aJ a)(a^ a)a 2 structure which can be rewritten as a 
sum of (at) 4 a 4 and (a'l') 3 a 3 . The latter contributes to the 
three-body sector, though naively the diagram appears to 
affect only the four-body sector. Note also that a con- 
traction of the lines amounts to an additional loop, thus 
the latter contribution actually comes from the two-loop 
diagram. 

When the diagram with the additional loop is con- 
sidered, we should treat it in a symmetric and Galilean 
invariant way. Because of symmetry between the iden- 
tical particles, we should treat the two internal lines of 
the additional loop symmetric. Furthermore, we need to 
impose the cutoff on the relative momentum of the two 
lines to maintain Galilean invariance. In this way, we 
realize that the "total momentum" of the loop must be 




FIG. 3: A time-ordered diagram which contributes to four- 
body local counterterms when the loop momentum is set 
in the shell. The resulting counterterms are not normal- 
ordered. The rewriting them in the normal-ordered form gen- 
erates normal-ordered three-body operators as well as normal- 
ordered four-body operators. 



A/2 + k 




FIG. 4: The diagram obtained from Fig. by normal- 
ordering. The external energies and momenta are set to zero 
for simplicity. The momenta are assigned so that Bose sym- 
metry and the Galilean invariance is maintained. Note that 
all the momenta of the lines are related to the shell-mode 
momentum. 



the shell mode, while cutoff is imposed on the relative 
momentum, as shown in Fig. [4j 

We go on to the third example in Fig. \b\ The diagram 
has six incoming and six outgoing external lines. Naively 
it only produces effective six-body operators. However, 
since it is not normal-ordered, rewriting it in the normal- 
ordered form is necessary to obtain the correct running 
of the coupling constants. There are several ways of con- 
tracting the lines. In particular, it gives rise to three- 
body operators, which actually come from four-loop di- 
agrams because the contractions amount to loops. In a 
similar way, we can continue to one-loop diagrams with 
an arbitrary number of external lines. Again, Bose sym- 
metry and Galilean invariance force a particular momen- 
tum assignment so that the momenta of the lines of the 
additional loops are related to the shell mode. 

There is another kind of diagrams generated by 
normal-ordering, in which the contractions give rise to 
internal lines which get across the shell-mode loop. Sec 
Fig. [6] for examples. Such diagrams do not contribute 
the RGEs. In the usual relativistic field theory, such 
diagrams can be obtained after constructing the effec- 
tive local vertex by connecting the legs. That is one of 
the reasons why the RGEs have the "one-loop" prop- 
erty mentioned above. In the nonrelativistic field theory, 
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FIG. 5: A time-ordered diagram which contributes to six- 
body local counterterms when the loop momentum is set in 
the shell. The resulting counterterms are not normal-ordered. 
The rewriting them in the normal-ordered form generates 
normal-ordered three-body operators, as well as operators 
which act only on other sectors. All the momenta of the 
lines of the additional loops are related to the shell mode. 




FIG. 6: Examples of the diagrams in which the additional 
lines get across the shell-mode loop. The momentum associ- 
ated with these lines are not related to the shell mode. 



because there is no antiparticle, tadpole-type diagrams 
cannot occur. We can just disregard these contractions. 
The crucial difference between the diagrams considered 
in the previous paragraphs and the ones with lines across 
the shell-mode loop is that in the former the additional 
lines are forced to related to the shell mode, while in the 
latter they are not. 

The analysis given above provides the solution of the 
"puzzle." The one-loop property persists even in the 
three-body sector, but the generated diagrams are not 
in general in the normal-ordered form. In order to ob- 
tain the running of the coupling constant for the oper- 
ators in the Lagrangian of a nonrelativistic theory, one 
needs to rewrite the contributions in the normal-ordered 
form because the operators in the Lagrangian are usually 
assumed to be normal-ordered. In this rewriting, some 
contractions may be needed. A contraction amounts to 
an additional loop. The diagrams obtained by normal- 
ordering in which all the lines of the additional loop are 
forced to be related to the shell mode, with an arbitrary 
number of loops, may contribute to the running of three- 
body operators. 

Note that no multi-loop contribution can occur in the 
two-body sector. This is a special feature of the two-body 
sector. On the other hand, it is much more involved 



in many-body sectors. In the three-body sector, as we 
will see, the infinitely many diagrams contributing to the 
running of the three-body coupling do can be rcsummed. 

IV. DERIVATION OF RGES 

As we explained in the previous section, there are 
multi-loop contributions to three-body operators ob- 
tained by normal-ordering rewriting of the one-loop dia- 
grams that are not in the normal-ordered form. It turns 
out that there are five types of such "normal-ordered one- 
loop" (NOOL) diagrams, together with two genuine (i.e., 
already normal-ordered) one-loop ones. They are given in 
Fig. [7J Note that the multi-loop contributions are neatly 
written in terms of the two-body scattering amplitude (a 
blob in Fig. [Jj which satisfies the Lippmann-Schwinger 
equation depicted in Fig. [5] Note also that the diagram 
in Fig. [5] appears as a part of the first diagram in Fig. [JJ 




^0^ )Q( 

FIG. 7: The NOOL diagrams together with the two genuine 
one-loop diagrams. The shaded blob stands for the two-body 
amplitude defined in Fig. [8] 




FIG. 8: The Lippmann-Schwinger equation for the two-body 
scattering amplitude, represented as a shaded blob. The ver- 
tex actually stands for a collection of vertices, cq, C2, and 62 
in the present approximation. 

The RGEs can be obtained by setting the loop mo- 
menta of the seven diagrams in the shell and canceling 
their cutoff dependence by adding suitable three-body 
counterterms, as we depicted in Fig. [9l 

It is important to notice that a sum of additional loops, 
as those in Fig. [SJ together with the two-body coun- 




FIG. 9: A diagrammatic representation for the RGEs of the 
three-body operators. Propagators with the momentum being 
in the shell are denoted by double lines. 
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terterms, forms the cutoff-indcpcndcnt two-body ampli- 
tude so that it does not contribute to the running of the 
coupling constants of the three-body operators. Simi- 
larly, the cutoff dependence of the additional loops on 
the oblique sides of the trapezoids in the second, third, 
and the fourth diagrams in Fig. [5] is already taken into 
account by the sixth and the seventh diagrams in Fig. [9] 
In this section, we consider only the terms that are 
explicitly shown in Eq. (|2.1[) . that is, two-body opera- 
tors up to including 0{p 2 ) and a three-body operator 
without derivatives. This truncation makes our RGEs 
approximate. Higher order corrections are discussed in 
Sec. EH 

The two-body amplitude can be treated separately. We 
assume that the cutoff A is small enough so that the two- 
body amplitude is given in the effective range expansion 
form. The renormalizcd off-shell amplitude may be writ- 
ten as 



— -/i+ir e /i 2 -— (2/i 2 + fc 2 + fc 2 
a 2 2 aa 



(4.1) 

where a 2 and r e are the scattering length and the effective 
range respectively, k\ and k 2 are incoming and outgoing 
relative momenta, 



/'• 



(4.2) 



with = (P°,P) being the total momentum. (Note 
that Eq. (|4.ip is different by a symmetric factor 1/2 
from that for the "spinless nucleon" given in Ref. [26[ 
for which, to mimic the nucleon case, we intentionally ig- 
nore the factor.) We have an off-shell parameter h com- 
ing from the redundant operator. The parameters a 2 , r e , 
and h are cutoff independent. 

The RGEs for the coupling constants for the two-body 
operators can be obtained in Ref. (26j : 



d -§ = (1-X)(Y + 3X 2 ) /X 



(4.3) 

2XY - Y) /X 2 , (4.4) 



— = -Y(6X 3 -5X 2 
dt v 

— = - (&X Z Z - 3X 2 Z + 2XYZ + Y 2 ) /X 2 ,(4.5) 
dt 



where X, Y, and Z are combinations, 



X = l + -(y + z), 



Y 



- 5 (y + *) 2 , 



Z = 2y+-(y + z) 2 , 
of the dimensionless coupling constants, 



MA 
4^2 



co, y 



A/A 3 

1^2 



■4c 2 



4tt 2 



bo. 



(4.6) 
(4.7) 
(4.8) 

(4.9) 



The parameter t is defined as t = ln(Ao/A), where Ao is 
the physical cutoff which is the limit of the applicability 



of the EFT, and A is the floating cutoff. There is a non- 
trivial fixed point, (X* ,Y* , Z*) = (1,-1,-1), besides 
the trivial one, (1, 0, 0). 

These RGEs have an analytic solution: 



X 
Y 
Z 



(CA)- 1 - 1 



(CA)- 1 - 1 -CA 2 ' 

(CA)' 1 - 1 
[(CA)- 1 - 1 - CA 2 ] 5 

C"A - 1 
[(CA)- 1 - 1 - CA 2 ] 5 



(4.10) 



where C, C, and C" arc the integration constants. They 
arc related to the effective range expansion parameters, 
a 2 , r e , and h, by 



C 



2a 2 

7T 



C = h, C" 



(4.11) 



In terms of them, the three-body RGEs are written as 
follows: 



dv 
~dl 



1 

s 2 



T-- b U\-c 



12b 



S 2 



T- 3 



U - 4 - 6V 



V 2 

+ \2—{T-W) 

where we have introduced several combinations, 
S = X 2 -aY, 



T 



YS + -{X 2 Z + Y 2 



U = X (X - 1) (X 5 



Y) 



V = Y + 3(X -1) + -(X -l) 2 , 
5 

together with numerical constants, 
a = — tt-1, 



(4.12) 

(4.13) 
(4.14) 
(4.15) 
(4.16) 



V3 

7 

4 

V3 



V3 . V3 

arctan — , 

2 2 ' 



7V3 ^ V3 

7r arctan — 

3 4 18 2 



V ( 4 - 17 ) 



Note that in these equations, X, Y, and Z are the 
solution, Eq. (|4.10[) . with the integration constants, 
Eq. gTLj. 

In Appendix [B] a sample calculation of a contribution 
to the RGE is given for the purpose of illustration. It 
is a bit difficult to find out the shell mode contribution 
of the two-loop diagram shown in the fifth diagram in 
Fig. [3] In Appendix [C] we explain how to identify the 
shell mode contribution. 

It is important to note that the RGE now has terms 
proportional to v 2 . The appearance of these terms allow 
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FIG. 10: The dimeron representation of the two-body ampli- 
tude, composed of the dressed dimeron propagator and the 
two-particle-dimeron vertices at the ends. 



d6 



ex 




the periodic behavior of v in the RG evolution, as we will 
see in the next section. 

It is interesting to note that these seven diagrams in 
Fig. [7] are the irreducible building blocks of which the 
three-body amplitude is composed, just as simple two- 
point bubbles are the building blocks of the two-body 
amplitude. 

The amplitude composed of these blocks can be com- 
pared with that obtained in the formulation with the 
dimeron field. The two-body amplitude is equal to the 
(dressed) dimeron propagator multiplied by the factor 
g\ coming from the two-particle-dimeron vertices at the 
ends. See Fig. [TU] With this identification, one easily 
sees the correspondence between the three-body ampli- 
tude in our formulation and that in the formulation with 
dimeron. 

There is however a difference: in our formulation, a 
one-loop diagram with a two-body vertex and a three- 
body vertex (which we call a "two-three loop") con- 
tributes, while in the existing calculations with dimeron 
these contributions are not included. 

Note that the cutoff dependence of the two-three loop 
cannot be renormalized as a part of two-body amplitude 
because the two-body amplitude is composed solely of 
two-body vertices, but is canceled by three-body coun- 
terterms, as we explained earlier. See the sixth and the 
seventh diagrams in Fig. |H1 The momentum dependence 
becomes important when the two-three loop is embed- 
ded in the three-body amplitude. When the shell-mode 
momentum is considered as in the second, third, and the 
fourth diagrams in Fig. |H1 the momentum dependence of 
the two-three loops becomes an additional source of the 
cutoff dependence. 

In the existing calculations with dimeron, the D<j)- 
D<j) vertex is momentum independent even for the full 
dimeron lines. On the other hand, as depicted in Fig.QTJ 
a single insertion of the three-body vertex in the two- 
body amplitude cannot be represented by a momentum- 
independent Dfjj-Dfjj vertex, but has momentum depen- 
dence due to the two-three loops. 



FIG. 11: The insertion of the three-body vertex in the two- 
body amplitude. In the third line, the two-body ampli- 
tudes are replaced by the full dimeron propagator. Because 
of the two-three loops, it is not equal to the (momentum- 
independent) D<f>-D<f) vertex. 



RG FLOWS AND THE EFIMOV 
PARAMETER 



A. Flows in the r e = h = subspace 

Let us begin, for simplicity, with the RG flows on which 
the effective range, r e , and the off-shell parameter, h, are 
zero. In this case, the RGEs are drastically simplified 



and are written in terms of £ = Y/X 



as 



dt 
d 

—i 

dt 



-£(! + £). 

12£ 3 / 126£ 2 



1 



36 2 g 
l-a£ 



at 



6t)v 



(5.1) 



(5.2) 



Note that the RGE for £ is the same as that for x in 
the leading order approximation. The RG flow in the £-v 
plane are shown in Fig. 1121 The line £ = — 1 corresponds 
to the nontrivial fixed point of the two-body RGEs. The 
flow there goes down periodically, exhibiting the limit 
cycle behavior. In addition to the trivial fixed point, 
we find a nontrivial fixed point (£*,«*) = (0, — 4/c) = 
(0, —8.126 ■ ■ • ), which we call the Borromean fixed point, 
though it apparently has nothing to do with the known 
Borromean systems. In the full set of RGEs, the Bor- 
romean fixed point is at (x*, y*, z*, V+) = (0, 0, 0, — 4/c), 
so that there are no two-body interactions but only the 
three-body interaction. 

One may notice that there is a singularity in Eq. (|5.2[) 
at £ = 1/a. This is due to the pole of the two-body 
amplitude which corresponds to the bound state. See 
Appendix |D] for more details. 
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FIG. 13: The running of v and £ as functions of t in the two- 
body strong-coupling phase (£ < —1). The initial value is 
taken as v(t = 0) = and £(t = 0) = -1.001. 



FIG. 12: The RG flow in the £-v plane. The arrows indicate 
the directions of the flow to the infrared (larger t) . 



B. Limit cycle 

The limit cycle behavior is due to the v 2 term in the 
RGEfor v. With (x,y,z) = (-1,-1/2,1/2) ((X,Y,Z) = 
(1, —1, —1)), where, = oo and r e = h = 0, the RGE 
for v becomes 



dv 
~dt 



= A + Bv + Cv 2 , 



with 
A 



48 



D 



486 



C 



lib 2 



so that the discriminant is negative, 



D = B — AAC 



-4.27374- •• < 0. 



Thus the solution is given by 



—B 



-(t + to) 



2A 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



where to is a constant of integration. It shows the peri- 
odic behavior with the period 



T 



2 71 



3.03932. 



(5.7) 



In terms of the Efimov parameter, Sq, it is written as 



s = -~ 1.03365, 
which is only 2.7% off the Efimov's value. 



(5.8) 



It is interesting to compare the value with that ob- 
tained in the "point-like" approximation in the dimeron 
formulation [11], sq — 1.393, which is about 40% off the 
Efimov's value. The difference is due to the contributions 
from the two-three loops, as we explained in the previous 
section. 



C. Off the critical case 

We can investigate the off-critical, i.e., finite scattering 
length, behavior. Fig. Q2] and Fig. [TJ] show the running 
of v and £ as functions of t starting with £(i = 0) + 
1 = -0.001 and £(i = 0) + 1 = 0.001, respectively. If 
the flow is close enough to the critical line (£ = —1), v 
diverges (to negative infinity) a finite number of times 
before going to the trivial fixed point (for £ > —1) or to 
the positive infinity (for £ < —1). We suspect that the 
occurrence of divergences corresponds to the existence 
of three-body bound states, as for the critical case in 
which the flow diverges infinite times, corresponding to 
an infinite number of bound states. Thus, if the flow is 
close enough to the critical line, there are a finite number 
of three-body bound states. 

The behavior in the two-body strong-coupling phase 
(£ < — 1) is particularly interesting from the NEFT point 
of view. In the 3 Si- 3 -Di channel, the two-nucleon sys- 
tem is in the strong-coupling phase and very close to 
the critical case. There is a shallow two-body bound 
state, deuteron. The binding energy of the triton (E t = 
8.48 MeV) which is in the spin one-half channel in the 
neutron-deuteron system, is larger than that of deuteron 
(Ed = 2.2246 MeV). Even though the present theory is 
bosonic, we expect that the same RG structure persists in 
the fermionic theory, and that the interplay between the 
three-body and the two-body bound states of the bosonic 
theory may explain that of the nucleon systems. 
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FIG. 14: The running of v and £ as functions of t in the two- 
body weak-coupling phase (£ > —1). The initial value is takes 
as v(t = 0) = and f (t = 0) = -0.999. 

The coupling £ diverges to negative infinity at a finite 
value of t (£2) in the two-body strong-coupling phase. It 
implies the existence of the two-body bound state and 
the value t% would correspond to the binding energy of 
it. Once £ diverges to negative infinity, the flow appears 
from the positive infinity of £, decreasing to the trivial 
fixed point, £ = 0. (The continuity may be easily seen by 
introducing a change of variable from £ to 6 = tan -1 £. 
The RGE (JET) becomes d9/dt = -sin6>(cos6> + sin 6). 
The transition from the negative infinity to the positive 
infinity of £ corresponds to passing — ir/2 of 9.) 

When v diverges to negative infinity in the two-body 
strong-coupling phase (£ < —1), it occurs before £ di- 
verges. In the positive £ region, on the other hand, we 
numerically find, by examining the flow, that all the flows 
which come from the two-body strong-coupling phase 
never diverge to negative infinity. The above observation 
implies that, if both the two-body and the three-body 
bound states exist, the binding energy of the two-body 
bound state is always smaller than the three-body bound 
state. 



D. Borromean fixed point 

To our best knowledge, the existence of the Borromean 
fixed point has never been noticed in the literature. The 
existence of it and of the critical surface on which the 
Borromean fixed point resides implies a new phase where 
the two-body interactions do not support the two-body 
bound states, but the strong three-body interaction gives 
rise to the three-body bound states. 

We linearize the RGE near the fixed point, (£, v) = 

*(2) :)(£)■<»> 
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FIG. 15: The effects of nonzero values of h and r e . The red 
line is the solution Eq. (|5.6[) . with r e = h — 0. The green line 
is the case with r e Ao = 1 and h — 0, where Ao is the physical 
cutoff. The blue line is the case with r e Ao = /iAq = 1. Three 
of them have asymptotically the same period. 



By diagonalizing it, we see that the scaling dimension 
of the relevant coupling is +4. In the next section, we 
investigate the Borromean fixed point once again with 
higher order corrections. 



E. Effects of nonzero values of r e and h 

Going back to the full set of RGEs, Eq. (|4TT2|) to- 
gether with Eqs. (|4.10[) . we can investigate the effects 
of nonzero values of h and r e . The results are shown in 
Fig. [15] It is interesting to note that they affect the flow 
of d , especially the "period," even though their effects 
damp as t grows. Physically speaking, it implies that the 
ratios of the binding energies of the three-body bound 
states varies and asymptotically becomes the universal 
value e~ 2 */ s ° for the bound states accumulate toward 
the threshold. It is surprising to see that the parame- 
ter, h, corresponding to the redundant operator has such 
observable effects, even though it does not contribute to 
the two-body observables. Furthermore, it does not even 
appear in the off-shell amplitude, Eq. (|4.1j) . in the critical 
case (l/a2 = 0)! 



VI. HIGHER ORDER CORRECTIONS 

In this section, we concentrate on the cases in which 
the two-body sector is on the nontrivial and the triv- 
ial fixed points, and consider the effects of higher order 
three-body operators. Note that all the irrelevant two- 
body operators do not contribute in these cases. 
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We include the following three-body operators, 
d 2 



^higher 



48 



24 



-h.C. 



V 2 



-h.c. 



.1) 



as 



MA e , A 6 1 

w = ,„ n „ a 2 , w = /0 _ 9 ^ 9 e 2 , s = it + -w. (6.2) 



(2tt 2 ) 2 



(2tt 2 ) 2 



A. On the nontrivial fixed point 



and obtain the RGEs for the dimensionlcss coupling con- 
stants, v, s, and w, the latter two of them are defined 



On the nontrivial fixed point of the two-body sector, 
(x,y,z) = (-1,-1/2, 1/2), RGEs are given by 



dv 

777 

ds 
~dl 
dw 
~dt 



48 

V37T 

64 
3x/3tt 
576 486 



486 

V37T 

246 8 

V3?r 9V3tt 



V3tt v/3 7T 

"1 4 

4 



io - 24 



3 _]_ 

V3n V 4 + 12 



,s 2 + 



3c 



V3?r 



V^TT 

■ 2 



7 

120 



1 



V3n 



326 
3V3tt 



12^7T 

126 2 



V3tt 



+ c ) vs 



x/3tt 



(96 + 4) 



s + 6 



86 
V3tt 



(6 + 3d) 



(96 2 + 26 + 3d) 



vs 



„ / 126 2 



26 
x/3tt 



(6.4) 



3 ,SW ' 



(6.5) 



where the constants 6, c are given in Eqs. (|4.17l) and d 
and e arc defined as 



V3 
V3 



2 

-3 + ir 



V3 , Vs 

— arctan — , 
3 2 ' 

V3 V3 

— arctan — . 
9 2 



(6.6) 



Note that RGEs for v and s do not depend onio, so that 
they can be solved without solving the RGE for w. In the 
following we mainly consider the flows in the v-s plane. 

These RGEs can be solved numerically. The t depen- 
dence of the coupling constants, v, s, together with the 
leading order one for v, is shown in Fig. 1161 Note that 
after an initial transient region, they exhibit periodic be- 
havior. Note also that they do not diverge. The RG flow 
in the v-s plane is shown in Fig. 1171 

We numerically read off the period of the limit cycle 
as 2.9810, which corresponds to the value of Efimov pa- 
rameter s ~ 1.0539. The value is about 4.73% off the 
Efimov's value. Compared with the leading order result, 
it deviates a bit more from the Efimov's value, though it 
is still very close. It might suggest that the convergence 
is not monotonous. 

It is interesting to see that the limit cycle behavior 
emerges very differently from the leading-order one dis- 
cussed in the previous section. In the leading order, since 
there is only one coupling, the periodicity is not possible 
unless it gets through the infinity. On the other hand, in 
the higher order, it is possible for the flows to have a pe- 
riodicity while staying finite. It explains why the drastic 



V, S 100 




FIG. 16: The t dependence of the three-body coupling con- 
stants, ii and s. The blue line stands for the leading order 
solution of v. 



change of the behavior of the RG flows occurs when the 
higher order contributions are taken into account. 

Another difference comes from the existence of the 
transient region. In the leading order, exact periodic- 
ity appears from the beginning. Actually the structure 
of the RG flow in the higher order is much richer. First 
of all, there is a nontrivial fixed point, numerically found 
to be at (v+, s*, w+) = (—3.14,10.05,-60.73). As shown 
in Fig.[T71 the flows from this nontrivial fixed point grad- 
ually approach to the limit cycle. On the other hand, the 
flows outside the limit cycle loop are grouped in the two 
categories: some flows are directly attracted by the limit 
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FIG. 17: The limit cycle behavior in the v-s plane. The blue 
loop is the limit cycle. Inside the loop there is a nontrivial 
fixed point (red bullet). 
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FIG. 19: The RG flow on the two-body trivial fixed point 
(£ = 0). We find another nontrivial fixed point in addition to 
the trivial and Borromean fixed points. 




FIG. 18: The RG flows in the v-s plane in a wide view. An- 
other fixed point is shown together with the limit cycle loop. 



cycle and gradually approach to it, and the others go to 
infinity (and come back from the opposite side) and then 
approach to the limit cycle. The existence of these two 
categories is related to another nontrivial fixed point at 
(-135.52, 403.62, -1069.65). See Fig. M 

The eigenvalues of the linearized RGEs at the fixed 
point located inside the limit-cycle loop are numerically 
obtained as 0.3288±2.3212i and 1.9817. The eigenvector 
of the real positive eigenvalue is in the w direction. It is 
known that the imaginary part of the complex eigenval- 



ues, 9, gives approximately the period of the limit cycle, 

2vr 



T 



Im(9' 



(6.7) 



With 9 = 2.3212, we find T = 2.7069, which corresponds 
to the value of Efimov parameter sq = 1.1606. 



On the trivial fixed point 



As explained in Sec. IV Dl we find a nontrivial fixed 
point of the three-body coupling do on the two-body triv- 
ial fixed point. In this section, we examine the effects of 
the higher order operators on it. 

The RGEs on the trivial fixed point of the two-body 
sector are given by 



dv 
~d~b 
ds 
~dt 
dw 
~dt 



a 2 1 

-4v — cv — —vs — tttttS , 



1_ 

120" 



= —6s — cvs — — s 2 
6 



—6w + 6ev 2 + -s 2 + 3cvs — 2cvw — -sw.(6. 



The RG flow in the v-s plane is shown in Fig. Q1J] We 
find the Borromean fixed point at the same value of v. In 
the v-s-w space, the fixed point is at (— 4/c, 0, — 48e/c 2 ) w 
(—8.13,0,-74.80). This strongly suggests that the exis- 
tence of the Borromean fixed point is not an artefact of 
the restriction of the set of operators. 

By linearizing the RGEs at the Borromean fixed point, 
we find the eigenvalues, v, and the corresponding (unnor- 
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FIG. 20: The RG flow in the v-s plane with £ = -0.25. 
The Borromean fixed point and the trivial fixed point move 
together and eventually fuse and disappear. 



malized) eigenvectors, (6v,5s,5w); 



4 : 




Note that there are two relevant operators. 

We also find other fixed points 
at (-18.45,18.48,-58.81) and at 

(-247.09,693.75,-917.51). The former is located 
at the center of the whirl, and the latter, which is not 
displayed in Fig. Q1J1 is related to the fixed point shown 
in Fig. [HI 

Actually Fig.lEiltransforms into Fig. OH (or Fig.[T7J) if {; 

is treated as a parameter changing from zero (the trivial 
fixed point) to — 1 (the nontrivial fixed point approached 
from the weak-coupling phase). See Fig. to Fig. [5U 
The Borromean fixed point and the fixed point at the 
origin move together to fuse and disappear, while the 
whirl becomes enclosed by the limit cycle. This is an 
example of bifurcation. 



FIG. 21: The RG flow in the v-s plane with £ = -0.5. There 
is no Borromean nor trivial fixed point. 
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FIG. 22: The RG flow in the v-s plane with £ = -0.75. A 
large limit cycle loop emerges. 



VII. SUMMARY 

In this paper, we performed the Wilsonian RG analysis 
of nonrelativistic three-body systems of identical bosons 
without introducing the dimcron field. We explained why 
the multi-loop diagrams contribute to the RGEs and ob- 
tained the RGEs in the restricted space of operators. The 
periodic behavior in the case of infinite two-body scat- 
tering length is obtained with the value of the Efimov 
parameter being very close to the Efimov's value. The 



difference between our result and that obtained by using 
the dimeron field is clarified. 

We found very rich structures of the RG flows: beside 
the periodic behavior responsible for the Efimov effect, 
we found the Borromean fixed point, which, to our best 
knowledge, has never been noticed in the literature. It 
implies the existence of the three-body strong-coupling 
phase. We also investigated the effects of finite values of 
the effective range and off-shell parameters, and found 
that they affects the ratios of the binding energies of the 
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three-body bound states. 

We also try to explain that the binding energy of triton 
is larger than that of deuteron on the basis of the study 
of the RG flows in our bosonic system. 

We then extended our set of three-body operators to 
0(p 2 ) and considered the case of the two-body fixed 
points. On the two-body nontrivial fixed point, the limit 
cycle of the three-body couplings is realized as a loop of a 
finite size, and two nontrivial fixed points are found. On 
the two-body trivial fixed point there is a nontrivial fixed 
point around which a whirling of the RG flow occurs, in 
addition to the Borromean and the trivial fixed points. 

The three-body strong-coupling phase, with the phase 
boundary on which the Borromean fixed point resides, 
has a clear physical picture: it corresponds to the sit- 
uation in which there are no (or very weak) two-body 
interactions, but a strong (short-range) attractive three- 
body interaction among the particles. We do not know 
if such a situation is possible in the real world. 
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Appendix A: Use of a dimeron 

In this appendix, we illustrate a way of introducing a 
dimeron field in the path integral formalism, and empha- 
size that it is difficult to obtain the relations between the 
coupling constants of the original theory and those in a 
theory with a dimeron. 

Let us start with the (Euclidean) partition function, 



Z= dtfdcpe- J dix 



(Al) 



where Ce is the Euclidean version of the Lagrangian 
given by Eq. (|2.1[) . We insert the identity, 



1 = J dDUD e -f di * ADtD , 



(A2) 



where A is assumed to be a positive constant for simplic- 
ity. It can also contain derivatives. In such a case, the 
resulting free-part of the dimeron Lagrangian contains 
derivatives. 

We then make a change of variables, 

D = D' (1 + acpU +■■■)+ 2 (6 + c$<t> + ••■), (A3) 

£>t = (D'Y (1 + O 0V +••■)+ (^) 2 ( 5 + c <^ + ■■•)' 

(A4) 

where a, b, and c are real constants. Here and hereafter 
the ellipses denote the terms of 0's and t/^'s with deriva- 
tives. Note that we assign D and D' the particle number 



—2 and +2 respectively, with the assignment of the par- 
ticle number — 1 to </>. We may include terms of higher 
powers in D, and D\ but such terms do not con- 
tribute to the three-body sector, and are thus ignored. 
Under the change of variables, we have 

AD^D = A{D'^D' + 6A ((£>') V+kc.) 

+ 2aA(L>') t J DW 
+ (c + ab) A ((£>') W + h. c .) 
i 2 a (4^a\ 2 _l 26cA (^^ 3 



b A 



(A5) 



This is not the whole story, however. One needs to 
consider the Jacobian contribution unless one uses di- 
mensional rcgularization, 

dD ] dD = dD' f dD'J(<f>, (A6) 
The Jacobian may be written as an action, 



J(0,0t) 



exp 



d X ^Jacobian 



(A7) 



up to a numerical constant, with C jacobian being ex- 
panded in terms of local terms, 



■'Jacobian 



5 M cfJ 4> + Aa 2 0) 2 + Ba 3 <j)) 



(A8) 

with 6M, A, and B are regularization dependent real 
constants. The calculation of the Jacobian must be con- 
sistent with the regularization of the other parts, and it 
is difficult to ensure the consistency. See Ref. [27| for an 
example of consistent calculation of a Jacobian factor in 
the lattice regularization. Here we do not try to calculate 
it, since it does not affect the following argument. 

The first term in Eq. (|A8[) may be absorbed in the defi- 
nition of the mass. With the other terms being included, 
the total Lagrangian becomes 

C' E = CE + AiD^D' + bA^D^^ + h.c.) 
+ 2aA{D')^D'(j) i (j) 
+ (c + ab) A ((£>') VV + h.c.) 

+ (b 2 A + Aa 2 ) (^cj)) 2 + (26cA + Ba 3 ) (0V) 3 
+ • • • . (A9) 

Thus, to eliminate the terms (^(j)) 2 , (^(j)) 3 , and 
(£>') W + h.c. in the whole La grangian, we choose a, 
b, and c to satisfy the following relations, 



= (c + ab)A, 

= -^ + b 2 A + Aa 2 , 

= -— +2bcA + Ba 3 
36 

By eliminating b and c we have 



(A10) 
(All) 

(A12) 



-lB-A]a 3 + ^a+^ = 0. (A13) 
2 J 4 72 v ' 



15 



There is at least one real solution for a. Once the so- 
lution for a is obtained, the solutions for b and c are 
obtained easily. We have shown that, even in the pres- 
ence of the contributions from the Jacobian, the terms 
(</>'!' </>) 2 , (^(j)) 3 , and (D')^^ 3 + h.c. can be eliminated 
but, as we claimed before, the relations between the cou- 
pling constants of the original theory and those in the 
theory with the dimeron are not simple at all. 



Appendix B: An example of NOOL diagram 
contributions to the RGE for do 



2tt 2 



c + 4c 2 A 2 
A 2 

< [-A 2 /M + ie] 
-iM 2 



'2 A 2 



M 



A z A 



c Q + 4c 2 A" + -^A 2 I A 



VgA A 2 A 2 * 

2 'T'T 



V3A A 2 A 2 \ JA 

— T'T r A 



.2tt 2 



4*-^ (x + y + z) A 



\/?>A A 2 A 2 * 



rfA. (B2) 




FIG. 23: An example of the NOOL diagram contributing to 
the RGE for do. 



In this appendix, we illustrate how the NOOL diagram 
contributions are evaluated by showing an example. 
The amplitude of the diagram shown in Fig. l23l is given 

by 



-ICQ 



d 4 k 

-zAU,-,- 



i4c 2 k 2 - i — k 2 
M 



k° 



k 2 
2M 



IC 



-k° 



k 2 

2M 



2tt 2 



x A 



—r dk[c + 4c 2 /c^ + -=-k 



V3k k 2 k 2 

TT'T'T 



M 

k 2 



-k 2 /M 



(Bl) 



where the external energies and momenta are set to zero. 
The first factor comes from the two two-body vertices. 
A stands for the two-body amplitude (|4.1[) . correspond- 
ing to the blob in Fig. [23l We have integrated over A; 
by picking up the pole at k° = —k 2 /2M + ie by the con- 
tour integration on the upper half plane. Note that other 
poles and the cut in the amplitude is on the lower half 
plane. The cutoff is introduced for the relative momen- 
tum, in order to preserve the Galilean invariance of the 
theory 0]. 

The contribution to the RGE comes from the shell- 
mode part of the integral, A — dA < k < A: 



We have the same contributions from the interchanges of 
the external momenta, so that we multiply it with 3x3, 
obtaining 



36i-^- (x + y + z) 2 A 



'v^A A 2 A 2 



dA. (B3) 



The two-body amplitude in terms of AT, Y , and Z is given 
by 



A(ji, klk 2 ) 
4tt 2 " 



MA 



it [i X 2 A i 2 X 2 Z 

1 — H h — h 1 

2 A Y A 2 \ Y 2 



k 2 + kj\ X(X-1) (Xl_ 
~Y 



A 2 



Y 



(B4) 



The combination 1 + X 2 /Y is written in terms of the 
two-body scattering length, a 2 , as 



X 2 
Y 



7T 1 

2^A : 



(B5) 



and the combination A" 2 Z/y 2 + 1 is written in terms of 
the two-body effective range, r e , as 



X 2 Z 7T , 

T^ + 1= 4 reA - 



(B6) 



Note that the cutoff A is small enough so that the effec- 
tive range expansion is valid. If we denote the two-body 

(2) 

cutoff A , above which the effective theory description 
breaks down, we assume 



A« A 



(2) 

■ 



(B7) 



We also assume that the effective range is of a natural 
size, i.e., 



r e A 2) ~ O(l), 



(B8) 



while the scattering length may (or may not) be fine- 
tuned. Thus, we have 



r e A < 1. 



(B9) 
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FIG. 24: The NOOL diagram, in which the identification of 
the shell mode contribution is not so obvious. 



The combination 3X(X - 1)(X 2 /Y + 1)/Y is related to 
the off-shell parameter h through 

X(X-l) fX 2 \ fU 2 , . 



We assume that h is also of a natural size, 

/>(a< 2) ) 2 ~o(i), 



thus 



hX 2 < 1. 



(Bll) 



(B12) 



Therefore the quantity in Eq. (|B10[) is much smaller than 
the one in Eq. (|B5|) . 

We may expand the corresponding terms in 
„4(^A/2,A 2 /4,A 2 /4), keeping 1 + X 2 /Y in the 
denominator. We finally get 



(2tt 2 ) 2 V 2 

where S, T, U, and V are defined in Eqs. (|4.16p . 



(B13) 



Appendix C: Shell mode contribution of the 
two-loop diagram with three symmetric lines 

In this Appendix, we consider the shell mode contribu- 
tion of the two-loop diagram shown in the fifth diagram 
in Fig. M 

This NOOL diagram is obtained by contracting two 
legs of the one-loop diagram as shown in Fig. [Ml The 
additional line connects the same vertices as two lines in 
the original diagram do. As we explained in Sec. IIIIB1 
the three lines should be treated in a symmetric way and 
we should impose the cutoff so that Galilean invariance 
is maintained. By doing so, the momenta of the three 
lines are related to the shell mode. 

A trouble is that there seems no obvious way to do so 
by just "modifying" the shell mode momentum assign- 
ment, as we did in Fig. 01 Fortunately, however, there is 
a trivial way to identify the shell mode contribution of 
the diagram starting from the amplitude, 

1(A) EE -(~id f J (^4/ (2f) 4_ fe 0_ fe2/2J/ + i£ 



k°/2 + l° - (k/2 + l) 2 /2M + ie 

i 

k°/2-P- (k/2-l) 2 /2M + ie 



(CI) 



where the external energy and momentum are set to zero 
for simplicity. The domain of the momentum integrations 
is restricted to the region |fc| < A and |Z| < A. These 
cutoffs do not break the Galilean invariance. Thus, the 
shell mode contribution is identified as 1(A) — I(A — dA). 
A simple calculation leads to 



1(A) - I(A - dA) 



6^Xcv 2 dA. (C2) 
A/A 5 v ' 



Appendix D: Pole at £ = 1/a 

In this appendix, we demonstrate that the pole at £ = 
1/a appearing in Eq. (|5.2[) is due to the existence of a 
two-body bound state. 

In Appendix |Bj we have shown how the two-body am- 
plitude is embedded in the three-body NOOL diagrams 
and contributes to the RGE for v. The two-body ampli- 
tude in the shell mode is given by 



(Dl) 



and we have argued that the effective range and the off- 
shell parameter terms can be expanded, so that effec- 
tively the amplitude appears as 




A 



A/A 



1 /i 
a 2 A A 



(D2) 



with [i = V3A/2. In the usual off-shell amplitude, [i is 
written as 



(D3) 



(D4) 



when the total momentum P is zero. On the other hand, 
because of Eq. (|B5|) . 



/i = ^-MP° + P 2 /A-ie, 
the value of \i corresponds to the energy 

P° = - 




1 _ 7T 1 

1 + e ~ 



£ = 1/a corresponds to 



(D5) 



(D6) 



(D7) 



which is nothing but the energy of a two-body bound 
state in the present approximation. 



1 


_V3 


a 2 A 


2 ' 


. flSJ 


, we get 


P° = 


1 




A/a 2 ' 
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1 *£<r 





























In the two-body RGEs, no value of the cutoff hits the 
pole of the amplitude, since /i is pure imaginary for any 
physical energy-momentum. In the present case, the two- 
body amplitude is embedded in the three-body diagrams, 
thus the "total energy-momentum" of the two-body am- 
plitude can take any real values so that /i can be real 
and hits the pole of the amplitude Eq. (jD2[) as the value 
CI2A changes. When £ comes close to 1/a, the three-body 
shell-mode contributions become huge, and the value of 
v changes rapidly to cancel the cutoff dependence. Near 
the £ = 1/a, the flows thus run almost vertically as shown 
in Fig. d5] 



FIG. 25: The RG flow in a wider region of the plane. A 
"wall" is seen at £ = 1/a. Rapid change of v near the "wall" 
is to cancel the enhancement due to the factor 1/(1 — <z£) in 
the RGE ([52J. 



There is however a "gate" , (1/a, — 2/ab), at which the 
right hand side of Eq. (|5.2[) is finite when £ approaches 
1/a. Apparently no flow can pass the £ = 1/a "wall" 
without getting through the "gate." 
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